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INTRODUCTION 


I . 


This  work  was  conducted  to  investigate  the  properties 
of  junction  switching  transistors  as  they  are  driven  into 
second  breakdown.  Second  breakdown  is  an  instability  that 
occurs  when  a  transistor  is  subjected  to  simultaneous  high 
collector  voltages  and  currents,  such  as  those  that  may  be 
induced  in  rhe  presence  of  an  electromagnetic  pulse.  An 
understanding  of  the  physical  mechanisms  which  lead  to 
second  breakdown  is  essential  for  the  prediction  of  circuit 
failure  in  such  an  environment.  Reverse  bias  second  break¬ 
down,  which  is  particularly  severe  because  of  the  high  power 
which  may  be  dissipated  in  the  device,  occurs  when  the  load 
is  inductive  and  is  initiated  when  a  transistor  is  being 
turned  off.  It  is  characterized  by  a  rapid  fall  in 
collector-emitter  voltage,  accompanied  by  an  increase  in 
collector  current,  and  generally  results  in  destruction  or 
degradation  of  the  device.  However,  a  nondestructive 
measurement  technique  has  been  developed  (Ref.  1)  in  which 
the  transistor  is  protected  by  way  of  a  collector  current 
crowbar  which  triggers  when  the  instability  has  been 
detected.  Such  a  system  was  constructed  and  used  at  Texas 
Tech  University  previously  (Ref.  2).  The  present  work, 
which  is  based  on  these  earlier  results,  had  a  two-fold 
objective:  to  investigate  in  more  detail  the  possible 
exclusion  of  a  thermal  initiation  mechanism,  and  to  extend 


It  had  two  subsidi- 


the  measurement  to  lower  temperatures, 
ary  objectives  as  well;  to  establish  the  limitations  of  the 
measurements  which  depend  on  the  test  system  itself,  and  to 
develop  a  system  for  the  rapid  acquisition  and  processing  of 
large  quantities  of  data. 


II.  EXPERIMENTAL  PROCEDURES 


GENERAL  DESCRIPTION  OF  THE  TESTER 

The  test  system  was  described  in  detail  in  Reference  2, 
but  a  review  of  its  operation  will  be  useful  here.  The 
transistor  under  test  (TUT)  is  turned  on  initially  with  a 
base  current  pulse  generated  by  the  forward  pulse  generator 
(Figure  1);  the  amplitude  of  this  pulse  is  controlled  by  the 
variable  gain  current  amplifier,  and  its  duration  by  the 
timing  section.  The  load  is  a  1  mH  inductor,  linear  to 
about  30  A.  The  collector  supply  is  a  variable  current 
limited  voltage  source  which  is  gated  on  and  off  by  the 
timing  circuit  of  the  base  drive.  The  initial  state  of  the 
TUT  is  controlled  by  varying  the  forward  pulse  amplitude, 
the  pulse  duration,  or  the  collector  supply  voltage. 
Immediately  following  the  forward  base  pulse,  the  timing 
circuit  triggers  the  reverse  pulse  generator,  which  drives  a 
second  variable  gain  current  amplifier.  After  the  charge 
stored  in  the  base  region  has  been  removed  by  the  reverse 
base  current,  the  collector  current  begins  to  fall;  the 
collector-emitter  voltage  increases  because  of  the  energy 
stored  in  the  magnetic  field  of  the  inductor.  At  some  point 
during  transistor  turn-off,  an  instability  generally  occurs, 
and  the  collector-emitter  voltage  falls  rapidly,  typically 
530  V  in  25  ns.  The  rapid  transition  causes  a  displacement 
current  in  the  detection  capacitor,  triggering  the  high 
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speed  crowbar.  The  collector  current  is  then  diverted 
through  the  crowbar  to  the  -140  V  power  supply.  Clamp 
supplies  limit  the  base  emitter  and  collector-emitter  vol¬ 
tages  (the  collector-emitter  clamp  is  not  shown). 

MEASUREMENT  APPARATUS 

Waveforms  were  recorded  on  a  Tektronix  7834  400  MHz 
storage  oscilloscope,  using  7A19  single-trace  600  MHz  verti¬ 
cal  amplifiers.  Two  time  bases,  the  7B80  and  7B85,  provided 
the  sweep.  The  collector-emitter  voltage  was  measured  with  a 
Tektronix  P6  057  1  00X  probe  having  a  1.4  GHz  bandwidth,  and 
the  base-emitter  voltage,  with  a  Tektronix  P6056,  10X, 

3.5  GHz  probe.  Base  and  collector  currents  were  measured 
with  Pearson  Model  411  current  transformers,  with  rise  times 
of  10  ns  and  outputs  of  50  mV  per  A. 

COLD  FINGER  AND  DEVICE  FIXTURE 

A  copper  cold  finger  was  constructed  to  permit  low 
temperature  cooling.  The  cold  finger  was  made  from  a  3/8 
inch  diameter  copper  rod  bent  into  the  form  of  a  J.  A  flat 
section  was  milled  on  the  short  end  and  a  hole  was  drilled 
to  permit  the  cold  finger  to  be  bolted  directly  to  the  TUT. 
A  small  circular  di3k,  1-1/2  inches  in  diameter,  was  welded 
onto  the  long  end  of  the  cold  finger;  this  increased 
surface  area  immersed  in  the  coolant,  in  this  case,  liquid 
nitrogen.  As  much  of  the  cold  finger  as  was  practical  was 
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previous  work  (Ref.  2),  although  still  in  a  nearly  loga¬ 
rithmic  sequence,  were  not  exactly  at  the  values  specified, 
particularly  at  high  currents.  (For  some  of  the  measure¬ 
ments  at  low  temperatures  in  this  work,  ,  the  reverse  base 
drive  could  not  supply  more  than  2  A  to  the  transistor.  It 
is  not  clear  what  caused  the  limitation.) 

TRANSITION  TIME 

The  transition  times  between  the  forward  and  reverse 
base  pulses  varied  with  the  pulse  amplitude,  from  around  100 
ns  at  equal  forward  and  reverse  base  drives  of  50  mA  to  some 
2  us  at  equal  base  drives  of  11  A;  Figure  9  illustrates  the 
latter  transition.  (Transition  time  measurements  were  made 
with  a  current  transformer  on  a  transistor  with  the  col¬ 
lector  open  circuited.)  This  is  an  important  limitation  of 
the  tester,  because  at  high  reverse  base  currents  the 
transistor  may  go  into  second  breakdown  before  the  reverse 
base  current  reaches  its  final  value.  Various  attempts  were 
made  to  improve  the  transition  time,  but  these  yielded  only 
minor  improvements.  This  problem  had  been  recognized 
previously  (Ref.  3),  but  it  had  been  assumed  that  the  turn¬ 
off  delay  time  was  long  enough  so  that  the  reverse  base 
drive  would  reach  its  final  value  before  second  breakdown. 
This  assumption  was  found  to  be  invalid  for  much  of  the 
present  work,  and  additional  work  will  require  redesign  of 
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EXPERIMENTAL  RESULTS 


TEST  CONDITIONS 

The  voltage  set  at  the  collector  supply  was  arbitrarily 
chosen  to  be  10  V.  The  emitter  base  voltage  was  clamped  ac 
6.5  V  and  the  collector-emitter  voltage  was  left  unclamped. 
The  forward  and  reverse  pulse  widths  were  set  to  4.5  and  10 
ms,  respectively.  The  former  is  long  enough  to  establish 
the  initial  state  of  the  transistor,  but  short  enough  to 
prevent  saturation  of  the  collector  current  transformer. 
Measurements  were  made  at  room  temperature,  approximately  20 
to  25  °C,  and  at  low  temperatures  between  -30  and  -60  °C. 
The  forward  base  current  pulse  amplitude  was  held  constant 
at  some  value,  while  the  reverse  base  current  pulse  ampli¬ 
tude  was  varied  in  the  approximately  logarithmic  sequence 
0.05,  0.1,  0.2,  0.5,  1,  2,  5  and  11  A.  At  first,  as  in  the 
previous  work  (Ref.  2),  the  current  levels  were  obtained  by 
setting  calibrated  controls  on  the  front  panel.  These 
controls  were  calibrated  by  inserting  a  1  ft  resistor  in 
the  base-emitter  terminals  on  the  test  fixture  and  measuring 
the  voltage  drop  across  the  resistor.  Later  experimen¬ 
tation  revealed  that  the  base  drive  was  load  dependent,  so 
the  calibrated  control  scheme  was  abandoned  in  favor  of 
using  a  current  transformer  with  a  fast  rise  time  and  lev 
percentage  droop  to  insure  chat  the  base  currents  were  r. 
t  h  e  specified  values.  Hence,  ohe  current  1 e  ;  a  1  ;  in  -  h - 
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Tha  value  fcr  instantaneous 
Vcq  at  the  final  paint  is 
663  V  at  6.  04  us. 

The  value  fcr  instantaneous 
Ic  ct  the  final  point  is 
,  731  A  at  6.  04  us. 

The  value  for  instantaneous 
POWER  ct  the  final  point  is 
534  W  at  6.  04  us. 

The  value  for  instcntaneous 
ENERGY  ct  the  final  point 
is  5370  uJ  at  6.  04  us. 


The  tima  at  10% 

of  the  second-breakdown  voltage 
is  To*.  622  us. 

The  value  of  instantaneous 
Vce  at  To  is  68.  3  V. 

Tha  value  of  instantaneous 
Ic  at  To  is  3.  42  A. 

The  value  of  instantaneous 
ENERGY 

at  To  is  47.3  uJ. 

The  change  in  ENERGY  from 
To  to  second-breakdown  is 
5320  uJ. 

To  to  seccn.d-fcreckdown  is 
5.  42  us. 


Figure  G.  Data  printout. 
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lector  current  (triangular 
waveform)  and  base-emitter 
voltage  (lowest  waveform) 
as  functions  of  time. 


1 


emitter  voltage,  collector  current,  instantaneous  collector 
power,  and  energy  at  the  final  point  was  outputted.  Finally, 
the  total  energy,  EC3,  dissipated  at  the  collector  between 
Tq  and  second  breakdown,  and  the  difference  between  Tq  and 
the  final  point,  that  is,  the  time  to  second  breakdown, 
teD,  were  calculated  and  printed.  Figures  6,  7  and  8 
illustrate  the  operation  of  the  programs.  Figure  6  is  a 
photograph  of  representative  waveforms.  Figure  7  illustrates 
the  various  curves  obtained  from  the  program,  and  Figure  8 
shows  the  data  printout  corresponding  to  the  curves.  The 
scale  factors  can  be  adjusted  to  change  the  dimensions  of 
the  plots;  also  a  plot  of  time  dependent  variables  against 
each  other  (collector  current  vjs.  collector  emitter  voltage) 


can  be  obtained. 


graphics  tablet  and  the  points  representing  the  exterior 
corners  of  the  photographic  coordinates  were  manually 
digitized  and  fed  to  the  computer.  The  degree  of  rotation, 
if  any,  of  the  coordinates  relative  to  the  graphics  tablet 
was  calculated,  then  used  to  rotate  all  subsequent  data  as 
required.  Remaining  data,  consisting  of  waveforms,  were 
digitized  manually.  The  rotated  data  were  stored  in  an 
internal  register  until  data  entry  was  complete  and  the 
program  was  instructed  to  store  the  data  on  floppy  disk.  The 
data  were  stored  in  an  x-y  format  for  later  calculations . 
A  plot  of  the  stored  data,  and  a  printout  of  important 
device  data  (manufacturer,  type,  fabrication),  device 
number,  test  conditions,  and  summary  results  were  provided. 

After  two  designated  waveforms  (relative  to  the  same 
time  scale)  had  been  entered,  additional  points  were  calcu¬ 
lated  by  interpolation,  and  a  product  curve  was  obtained. 
Here,  collector  current  and  collector-emitter  voltage  were 
entered  vs.  time,  and  instantaneous  power  was  calculated. 
The  time  integral  of  the  power  curve  was  then  obtained  to 
provide  the  energy  dissipated  at  the  collector  as  a  function 
of  time.  The  calculation  also  provides  the  time,  T  ,  at 
which  the  principal  input  curve  (collector-emitter  voltage) 
is  at  10%  of  its  peak  value;  collector-emitter  voltage, 
collector  current,  and  energy  at  this  time  were  then  out¬ 
putted.  Time  at  second  breakdown  (discontinuous  fail  in 
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Data  entry  and  analysis  system. 


Ill . 


DATA  ENTRY  AND  ANALYSIS 


The  planned  system  for  acquiring  fast  transient  data 
rapidly  from  large  numbers  of  devices,  and  for  storage  and 
analysis  of  these  data,  is  shown  in  Figure  5.  The  system 
controller  is  an  HP-85  desktop  computer  which  interfaces 
with  an  HP-2105  microcomputer.  The  microcomputer  buffers  a 
graphics  terminal,  a  printer-plotter,  and  a  magnetic  tape 
mass  storage  unit.  The  HP-85  itself  connects  to  a  multi¬ 
plexer  and  digital  voltmeter,  a  digital  oscilloscope,  a 
transient  digitizer,  graphics  tablet,  disc  storage  unit,  and 
printer-plotter.  The  three  instrument  groups  (HP-3497A  and 
3437A,  Nicolet  2090-III,  and  Tektronix  7612D)  represent  a 
capability  for  sampling  rates  varying  between  several  hours 
(HP-3497A  and  3437A)  and  5  ns  (Tektronix  7612D).  Because  of 
its  complexity,  it  was  not  possible  to  complete  the  assembly 
of  the  entire  system  during  the  contract  period;  existing 
and  projected  links  are  shown  in  Figure  5. 

It  was  not  possible  to  deliver  data  directly  from  the 
test  system  to  the  HP-85  via  the  Tektronix  7612D.  However, 
although  Intermediate  data  storage  on  film  was  required, 
parts  of  the  system  could  be  used  for  entry  and  analysis  of 
large  quantities  of  data.  Two  separate  processes  were 
required;  data  entry  using  the  program  MASTER,  and  calcula¬ 
tions  with  the  program  MATHPK.  (These  programs  are  listed 
in  Appendixes  A  and  B.)  A  photograph  was  affixed  to  the 


only  a  limited  number  of  devices  were  used.  Since  most  of 
the  previous  data  (Ref.  2)  were  obtained  with  Unitrode 
prototype  devices  similar  to  the  BUX  98  (30  A,  400  V), 
samples  of  these  transistors  were  also  used  here,  two  from 
one  lot  (NE)  and  two  from  a  second  lot  (NH).  Devices  19  and 
36  are  from  lot  NE-18W,  and  devices  26  and  37  are  from  lot 
NH-220-5  . 

The  difference  in  the  two  lots  is  the  emitter  design; 
Figure  3  is  a  photograph  of  the  NE  structure,  and  Figure  4 
shows  the  NH  pattern.  The  devices  are  double-diffused 
single  die  structures,  with  an  n+-epitaxial  substrate.  Base 
sheet  resistance  is  around  700  Q/a  under  the  emitter  stripe 
and  around  60  fi/n  between  emitter  fingers.  The  NE  die  has 
an  area  of  0.32  cm2,  the  NH  die,  0.22  cm2. 


(b)  Rear  view 


Figure  2.  Cold  finger  anti  device  fix¬ 
ture  . 
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wrapped  in  insulated  foam  and  duct  tape  to  minimize  thermal 
losses.  The  device  fixture,  a  section  of  printed  circuit 
board  for  the  minimization  of  strays,  was  coated  with  a 
silicone  encapsulant  and  a  sprayed  foam  insulator  to  prevent 
thermal  losses  and  inhibit  ice  buildup.  The  cold  finger  and 
device  fixture  are  shown  in  Figure  2.  However,  even  with 
insulation,  it  was  not  possible  to  reduce  the  temperature  of 
the  test  device  to  below  around  -60  “C.  There  was  not 
sufficient  time  during  the  course  of  this  work  to  design  and 
construct  a  low  stray  inductance  and  capacitance  fixture 
which  would  permit  measurements  at  lower  temperatures. 

A  Type  K  thermocouple  was  constructed  from  20  gauge 
chromel  and  alumel  wires;  the  bead  was  inserted  into  a  small 
hole  drilled  into  the  base  of  the  transistor  package  adja¬ 
cent  to  the  die.  The  thermocouple  was  held  in  place  with 
Epo-Tek  H70E,  a  thermally  conductive  and  electrically 
insulating  epoxy  adhesive.  The  thermocouple  leads  were 
connected  to  an  Omega  miniature  cold  junction  for  Type  K 
thermocouples.  The  output  was  measured  with  a  Keithley 
Model  177  digital  multimeter.  The  measured  voltages  were 
converted  to  temperature  using  standard  tables;  accuracy  is 
estimated  at  +  1/2  °C. 

TEST  DEVICES 

Because  of  the  extensive  measurements  which  were  per¬ 
formed  to  assure  characterization  over  all  operating  ranges. 


TYPES  OF  SECOND  BREAKDOWN 


In  general,  the  same  types  of 
in  the  previous  work  (Ref.  2)  were 
were  designated  as  Types  A,  B  and 
repeated  here  for  convenience. 


second  breakdown  observed 
also  observed  here.  These 
C;  their  descriptions  are 


Type  A  (Figure  10)  --  The  transistor  was  initially  in 
the  normal  active  region  or  just  quasisaturated.  This  type 
of  breakdown  was  normally  observed  when  both  forward  and 
reverse  base  pulse  amplitudes  were  small.  Following  a  delay 
time,  the  collector  emitter  voltage  rapidly  increased  to  the 


sustaining  voltage,  vcex(SUS)'  and  remained  there  for 


several  microseconds  before  increasing  to  the  second  break¬ 
down  voltage.  The  base-emitter  voltage  remained  positive 
during  the  delay  time  and  the  sustaining  period,  but  it 


decreased  as  VCE  increased,  approaching  zero  as  vc£  con¬ 


tinued  to  rise.  For  Type  A  breakdown,  t_Q  was  lonq  enough 

b  D 


to  allow  the  reverse  base  current  to  reach  its  final  value. 
The  collector  current  decreased  linearly  to  a  value  much 
less  than  its  maximum  at  turn-off.  Large  amounts  of  energy, 
of  the  order  of  several  thousand  microjoules,  were  delivered 
to  the  collector  during  Type  A  breakdown. 


Type  B  (Figure  11)  --  The  transistor  was  normally  in 
uasisatur union  trior  to  turn-off.  reached  v._..  ..... 


and  renal r. e d 


s  e  c  o  n  < 


reakdcwn.  The  t  - . 


Figure  11.  Typical  Type  E  waveforms. 

The  top  curve  is  1^,  (2  A 

per  large  division) ,  the 
center  curve  is  VRER  (0.5  V 

per  large  division) ,  and 
the  lowest  curve  is  V_„_. 

(100  V  per  large  division) . 
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emitter  voltage  remained  positive  prior  to  second  breakdown, 
and  the  base  current  generally  had  time  to  reach  its  final 
value.  The  collector  current  did  not  decrease  appreciably 
during  turnoff.  The  energy  delivered  to  the  collector 
during  Type  B  breakdown  was  also  on  the  order  of  a  few 
thousand  microjoules.  Type  B  was  found  to  be  relatively 
uncommon . 

Type  C  (Figure  12)  --  Type  C  breakdown  was  observed  for 
all  high  forward  base  currents,  and  for  low  forward  base 
drives  with  high  reverse  base  currents.  The  collector- 
emitter  voltage  did  not  reach  its  sustaining  value  and  the 
collector  current  did  not  decrease  appreciably.  The  base- 
emitter  voltage  became  negative  before  second  breakdown.  The 
energy  delivered  to  the  collector  during  Type  C  second 
breakdown  varied  from  a  few  thousand  microjoules  to  less 
than  a  hundred;  however,  the  instantaneous  power  at  the 
collector  at  second  breakdown  was  generally  larger  than  for 
Types  A  or  B. 

To  test  reproducibility,  the  measurements  made  in  the 
previous  work  (Ref.  2)  were  repeated.  A  minor  change  in  the 
base  drive  circuit  (Figure  13)  made  it  possible  to  increase 
the  forward  base  drive,  so  that  the  device  could  be  driven 
harder  toward  initial  saturation.  Only  Type  C  breakdown  was 
observed  with  the  higher  forward  base  drive  regardless  of 


Typical  Type  C  waveforms. 
The  top  curve  is  lc  (2  A 

per  large  division) ,  the 
center  curve  is  V  (0.5  V 

dLK 

per  large  division) ,  and 
the  lowest  curve  is 

(100  V  per  large  division). 


fication  for  higher  forward  base  drive 


the  amplitude  of  the  reverse  base  current  pulse. 

LOW  TEMPERATURE  MEASUREMENTS 

It  is  interesting  to  compare  the  low  temperature 
results  with  those  made  at  room  temperature.  In  general 
(although  not  always),  the  collector-emitter  voltages  at 
second  breakdown,  VSB<  were  lower  than  those  obtained  at 
room  temperature.  Both  the  delay  time,  t^,  and  the  time  to 
second  breakdown,  t^g,  decreased  at  the  lower  temperature, 
and  breakdown  shifted  from  the  slower  types  to  the  more 
rapid  types;  i.e.,  under  the  same  base  drive  conditions,  a 
device  which  exhibited  Type  A  breakdown  at  room  temperature 
displayed  Type  B  or  C  at  a  low  temperature,  and  a  device 
with  Type  B  at  room  temperature  shifted  to  Type  C  at  a  low 
temperature.  Table  1  contains  representative  data  comparing 
the  behavior  of  a  device  (Number  19)  at  room  temperature  and 
at  low  temperature.  Also,  the  devices  appeared  to  be  more 
vulnerable  at  lower  temperatures,  degrading  considerably 
with  repeated  measurement.  This  degradation  was  manifested 
as  a  rapid  increase  in  the  collector-emitter  voltage  during 
turnoff,  sometimes  as  great  as  400  V  in  10  ns.  This  was  not 
a  reproducible  phenomenon.  As  the  test  continued,  one  of 
the  devices  would  not  sustain  a  large  collector-emitter 
voltage  during  turnoff,  and  vn„  rose  to  only  10  to  15  v.  A 
faster  crowbar  or  reduced  clamp  voltages  may  be  required  to 
study  second  breakdown  at  low  temperatures. 
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Table  1 


COMPARISON  OF  ROOM  TEMPERATURE  AND  LOW  TEMPERATURE 
BEHAVIOR  FOR  DEVICE  NUMBER  19 


Note 


(forward  base  drive),  I0R  (reverse  base  drive). 


and  I___  (collector  current  at  second  breakdown)  are  in 

L  b  d 

amperes;  VeD  ( col lector -emitter  voltaqe  at  second  breakdown) 

O  D 


is  in  volts;  and  t. 


(time  to  second  breakdown)  is  in  micro¬ 


seconds  . 
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ENERGY  AND  SECOND  BREAKDOWN 


Voltage  and  current  waveforms  were  digitized  on  the 
graphics  tablet  and  stored  on  disks.  These  data  were  then 
processed  to  obtain  a  voltage-current  product,  which  was 
integrated  numerically  to  determine  the  energy  delivered  by 
the  load  inductance  to  the  collector  during  turn-off.  The 
energy  was  then  examined  as  a  function  of  the  time  to  second 
breakdown,  tCD,  which  is  defined  as  the  interval  between  the 

O  D 

time  at  which  the  collector-emitter  voltage  is  at  10%  of  the 
voltage  at  second  breakdown,  V_n,  and  the  time  at  which  V’ 
falls  abruptly.  The  delay  time,  t^,  is  the  interval  between 
the  time  at  which  the  monostable  multivibrator  generating 
the  reverse  base  current  pulse  is  triggered  and  the  time  at 
which  the  collector-emitter  voltage  increases  to  10%  of  its 
peak  value;  this  is  somewhat  different  from  the  usual 
definition  of  delay  time,  but  is  more  consistent  with  the 
time  to  second  breakdown. 

The  energy  vs.  te_  data  were  plotted  several  ways: 

O  D 

using  all  breakdown  data,  by  second  breakdown  type,  and  with 
forward  base  current  as  a  parameter.  Figures  14  through  25 
are  plots  of  energy  deposited  at  the  collector  at  room 
temperature  v_s.  time  to  second  breakdown  for  Device  Number- 
19.  These  are  quite  typical  of  the  data  for  the  other  three 
devices.  It  is  to  be  noted  that  Type  c  behavior  predomi¬ 
nates,  and  that  Type  B  is  uncommon.  As  forward  base  current 
is  increased,  the  energy  curves  are  shifted  up  somewhat. 
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Table  3. 

SUMMARY  OF  THE  RESULTS  OF  BEST  FITS  TO  THE  ENERGY  VS. 
TIME  TO  SECOND  3RSAKD0WN  DATA  FOR  DEVICE  NUMBER  26 


Plot 

a 

b 

All  Types 

'i  i  n  ; 

1 .00 

Type  A 

911 

1.32 

Type  B 

- 

- 

Type  C 

2910 

1.17 

IBp  =  0.05  A 

661 

0.95 

I BF  =  0-10  A 

1389 

0.92 

I BF  =  0-20  A 

2004 

0.97 

I gF  =  0-50  A 

2059 

1 .02 

^  BF  -  ^  ®  ^ 

3135 

1 .24 

I BF  =  2.0  A 

2987 

1.14 

1  BF  5-°  A 

3212 

1.02 

I BF  =  11*0  A 

34  52 

1.04 

-51- 


Table  2 


SUMMARY  OF  THE  RESULTS  OF  BEST  FITS  TO  THE  ENERGY  VS. 
TIME  TO  SECOND  BREAKDOWN  DATA  FOR  DEVICE  NUMBER  IS 


Plot 

a 

b 

Ail  Types 

2824 

1 .03 

Type  A 

1547 

1.09 

Type  B 

3431 

1.25 

Type  C 

3651 

1.12 

IBF  =  0-05  A 

1187 

0.98 

IBF  =  0.10  A 

2839 

1.03 

ID„  =  0.20  A 

BF 

2733 

1.01 

IBF  =  0.50  A 

2791 

1.02 

1 BF  =  1‘°  A 

4228 

1.12 

1 BF  =  2*°  A 

4055 

1.09 

I BF  =  5-°  A 

4  321 

1.08 

IBp  =  11.0  A 

4619 

1.08 
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This  effect  will  be  discussed  below.  The  highest  energies 
are  delivered  to  che  collector  for  Type  A  breakdowns;  Type  B 
breakdowns  represent  an  intermediate  energy  range,  and  Type 
C  breakdowns,  somewhat  lower  energies.  Because  device 
characteristics  changed  significantly  during  the  low  temp¬ 
erature  measurements,  it  was  not  practicable  to  digitize  and 
process  all  the  low  temperature  data.  However,  one  of  the 
devices  (Number  19)  appeared  to  have  undergone  less  degrada¬ 
tion  than  the  others,  so  a  partial  characterization  was  made 
for  it.  The  results  are  shown  in  Figure  26. 

The  results  of  the  best  fits,  Een  =  at^  ,  to  the  room 

bn  bo 

temperature  data,  are  summarized  in  Tables  2  through  5.  The 

best  fit  to  the  low  temperature  data  of  Figure  26  is  the 

relationship  E„n  =  34  51t^*  0  5  3  .  It  is  clear  that  the 

relationship  between  Esg  and  ts0  is  essentially  linear  for 

1/2 

all  cases.  Heating,  which  requires  a  t  dependence,  is 
excluded  as  an  initiating  mechanism  for  second  breakdown, 
at  least  for  the  Unitrode  devices,  and  perhaps  for  double- 
diffused  epitaxial  devices  generally. 

The  upward  shifts  in  the  energy  curves  which  were  noted 
above  were  examined  in  terms  of  the  dependence  on  forward 
base  current  of  the  coefficient,  a,  in  the  Eon  vs.  t_n  fits. 

b  d  b  d 

Figure  27  is  a  plot  of  a  vs.  I__;  these  data  were  fit  with  a 

r  b 

straight  line,  a  power  curve,  an  exponential,  and  a  loga¬ 
rithmic  curve.  The  data  were  best  fit  with  the  logarithmic 
curve,  a  =  a'  +  b '  1  r.  I  a  ;  the  results  of  the  fit  are  ore, - 
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Table  4. 

SUMMARY  OF  THE  RESULTS  OF  BEST  FITS  TO  THE  ENERGY  VS. 
TIME  TO  SECOND  BREAKDOWN  DATA  FOR  DEVICE  NUMBER  36  ' 


Plot 

a 

b 

All  Types 

2248 

1.00 

Type  A 

1168 

1.08 

Type  B 

- 

- 

Type  C 

3959 

1.28 

Ibf  =  0.05  A 

628 

1.03 

1 Bp  =  0.10  A 

1511 

1.01 

Igp  -  0.20  A 

2237 

1.05 

IfiF  =  °-5°  A 

2398 

1.04 

IBF  ~  ^ ’ 0  A 

4015 

1.18 

IBF  =  2-°  A 

3425 

1.01 

Ibf  =  5.0  A 

3468 

0.93 

Irp  =  A 

4  203 

0.96 
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Table  5. 

SUMMARY  OF  THE  RESULTS  OF  BEST  FITS  TO  THE  ENERGY  VS. 
TIME  TO  SECOND  BREAKDOWN  DATA  FOR  DEVICE  NUMBER  37 


Plot 

a 

b 

All  Types 

2387 

1.04 

Type  A 

1370 

1.10 

Type  B 

- 

- 

Type  C 

3014 

1.16 

Igp  =  0.05  A 

965 

1.09 

IBF  =  0.10  A 

1351 

1.03 

Ibf  =  0.20  A 

1764 

1.05 

IBF  =  0-50  A 

2213 

1.00 

:BF  =  1,0  A 

3374 

1.18 

:BF  =  2*°  A 

2667 

1.05 

1 BF  =  5*°  A 

3091 

0.89 

IBp  =  11.0  A 

34  23 

0.90 
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marised  for  the  four  devices  in  Table  6.  The  logarithmic 
fit  indicates  only  a  weak  dependence  of  the  energy  on 
forward  current;  what  is  more  interesting  is  the  relation¬ 
ship  between  a',  b'  and  the  mask  types.  The  respective 
constants  are  somewhat  larger  for  the  NE  masks  than  the  NH 
masks,  suggesting  that  a’  and  b'  may  be  related  to  emitter 
stripe  geometry  and  die  area. 


Table  6. 

SUMMARY  OF  THE  RESULTS  OF  LOGARITHMIC  CURVE 

FITS  TO  THE  COEFFICIENT,  a,  VS.  IQ„ 

or 


Device 

Mask  Type 

a' 

kj 

019 

NE 

3550 

562 

026 

NH 

2543 

499 

036 

NE 

2922 

639 

037 

NH 

2522 

458 

-56- 


V.  CONCLUSIONS 


The  room  temperature  results  indicate  that,  for  the 
devices  tested  here,  and  perhaps  for  all  double -diffused 
epitaxial  power  transistors,  a  thermal  initiation  mechanism  • 

for  second  breakdown  is  excluded.  This  is  an  important 
result,  and  its  generality  should  be  examined  by  testing 
similar  structures,  both  double-diffused  epitaxial  and  g 

triple-diffused,  from  other  manufacturers.  Comparisons 
should  be  made  with  single -diffused  transistors  (for  which 
Type  C  second  breakdown  should  be  confirmed)  to  verify  the 
linear  relationship  between  energy  and  time  to  second 
breakdown  at  these  temperatures,  and  particularly  to  study 
the  apparently  rapid  degradation  at  lower  temperatures. 

Finally,  the  data  entry  and  analysis  system  should  be 
completed,  to  permit  the  rapid  conditioning  of  the  large 
quantity  of  data  required  for  a  general  device  analysis,  and 
to  establish  procedures  and  criteria  for  acceptance  of  low 
vulnerability  devices. 
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APPENDIX  A 


Program  MASTER 


10  REM  * 

20  REM  * 

30  REM  *  MASTER  CONTROL 
40  REM  *  INTERFACE  ROUTINE 
50  REM  * 

80  REM  -Copyrights  9/26/83 
70  REM  -  gandalf  software,  inc. 

80  REM  *  Chuck  Graves.  Wizard 
90  REM  * 

100  REM  * 

10S  OPTION  EASE  19  CLEAR 
110  OUTPUT  706  ;  "IN" 

115  GOSUB  S6C0 

120  DISP  "How  many  curves  will  be  plotted?" 

130  INPUT  0 

140  IF  Q>0  AND  Q<7  THEN  170 

150  OISP  "Please  choose  a  number  1  to  6.  " 

160  GOTO  130 

170  INTEGER  T  <3>  .  P  <10>  .  U  <4>  .  V  <4>  .  X  (6.  1505  .  Y  <6.  150) 

175  DIM  P$ C2S3.  T* C503 

180  CLEAR  8  OUTPUT  706  i  "BP42.  1 25.  5" 

185  GOSUB  5300 

190  DISP  "Please  digitize  the  corners  of" 

200  DISP  "the  graph  grid  in  the  following" 

210  DISP  " orders  upper  left,  upper  right,  " 

220  DISP  "lower  left,  and  lower  right." 

230  FOR  1 1-1  TO  4 
240  GOSUB  1000 
250  IF  NOT  S  THEN  240 
260  U<I1)-T(1) 

270  V(I1)-T(2) 

280  NEXT  II 
295  GOSUB  4500 
300  FOR  J-l  TO  Q 

310  DISP  "Begin  digitizing  curve"i J» " 

320  Pl-J 

330  FOR  1-1  TO  150 
340  GOSUB  1900 

350  IF  PCI) -151  THEN  DISP  "Storage  file  is  full  <150  pts 

. "  9  WAIT  2000 

360  NEXT  I 

370  NEXT  J 

380  CLEAR 

390  OISP  "The  digitization  process  is” 

395  DISP  "complete.  If  you  would  like  to" 

400  DISP  "do  anything  else,  please" 

410  OISP  "digitize  the  desired  softkey." 

420  GOSUB  1900 

430  IF  K-5  AND  H*-"PL0T"  THEN  460 

440  IF  K— 4  AND  H*-"PL0T”  THEN  460 

450  GOTO  390 
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460  CLEAR 

470  DISP  “Raload  the  plottar  with  paper.  " 

480  DISP  "Press  CONTINUE  whan  ready.  " 

490  PAUSE 

495  PLOTTER  IS  705 
500  CLEAR 

505  SCALE  0.12032.0.8831 

510  DISP  "Do  you  want  a  grid  an  tha  plot?"  9  DISP  " (Y/N) 

it 

520  INPUT  G$ 

530  DISP  "How  many  horizontal  divisions?" 

540  INPUT  U1 

550  DISP  3  DISP  "How  much  is  each  division?" 

560  DISP  "  (magnitude,  unit)  5  mV  is  5.  mV.  " 

570  INPUT  U2.  U2$ 

580  CLEAR 

590  DISP  "How  many  vertical  divisions?" 

600  INPUT  VI 

610  DISP  3  DISP  "How  much  is  each  division?" 

620  DISP  "(magnitude,  unit)  5  mV  is  5.  mV.  " 

630  INPUT  V2.  V2* 

640  CLEAR 

650  DISP  "What  is  the  horizontal  label?" 

660  INPUT  U$ 

670  DISP  3  OISP  "What  is  the  vertical  label?” 

680  INPUT  V$ 

685  LORG  6  3  CSIZE  (H* (U (2) -U (1) ) /8831*100-10)  /20.  .  5 
690  LOCATE  10.  H*  (U  (2)  -U  (1)  )  /8831*100-10.  10.  H*  (V  (1)  -V  (3)  ) 
/8831*100-10 
692  FRAME 

695  SCALE  0.  U1*U2.  0.  V1*V2 
700  IF  G*-"Y"  THEN  730 
710  AXES  U2.  V2.  0.  0.  1.  1.  5 
720  GOTO  731 

730  GRID  U2.  V2.  0.  0.  1.  1.  1 

731  FOR  1 1-1  TO  U1 

732  MOVE  1 1*U2.  -  (V2/10) 

733  LABEL  VAL*(I1*U2) 

734  NEXT  II 

735  LORG  8 

736  FOR  12-1  TO  VI 

737  MOVE  -  (U2/10)  .  I2*V2 

738  LABEL  VAL*(I2*V2) 

739  NEXT  12 

740  IF  Q<-1  THEN  758 

741  OISP  9  OISP  "What  is  tho  second  vortical"  9  DISP  "la 
be  1  ?" 

742  INPUT  VI $9  CLEAR 

743  DISP  9  DISP  "How  much  is  each  division" 

744  OISP  "on  the  second  vertical  axis?" 

745  OISP  "  (magntude.  unit)  5  mV  is  5.  mV.  " 
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746  INPUT  V3.  V3* 

747  LORG  2 

748  FOR  1 1-1  TO  VI 

749  MOVE  U1*U2-U2/10. I1*V2 

750  LABEL  VAL$CI1*V3) 

751  NEXT  II 

752  OEG  9  LOIR  90 

753  LORG  6  9  CSIZE  CH* CU <25 -U Cl) ) /6831*100-10) /14.  .  5 

754  MOVE  U1*U2*1.  1.  .  5*V1*V2 

755  LABEL  Vl$,  "  C“:  V3*i  ")  " 

756  RAD  a  LOIR  0 

758  LORG  6  9  CSIZE  CH*  CU  C2) -U  Cl)  ) /8831*1Q0-*10) /14.  .  5 

759  MOVE  .  5*U1*U2.  -CV1*V2/10) 

760  LABEL  U$:  "  C"j  U2$i  ")  M 
770  OEG  9  LOIR  90  8  LORG  4 
780  MOVE  -CU1*U2/10).  -  5*V1*V2 
790  LABEL  V$,  ”  C":  V2$:  ")  “ 

795  RAD  9  LOIR  0 

800  N-U1*U2/CUC2)-UC1)) 

810  0-Vl*V2/CVCl)-VC3)) 

820  FOR  1 1-1  TO  Q 

830  MOVE  CXCI1.  1)-UC1))*N.  CYCI1.  1)-VC3))*0 
840  FOR  12-1  TO  PCI1)-1 

850  DRAW  CXCIl.  I2)-U<1))*N.  CYCIl.  I2)-VC3))*0 
860  NEXT  12 

870  NEXT  II 

871  DISP  a  OISP  "Haw  many  linos  doas  tha  main" 

872  OISP  "titlo  havQ?" 

873  INPUT  T1 

874  LORG  4 

875  FOR  1 1-1  TO  T1 

876  DISP  8  DISP  "What  is  lino"; II:"?"  9  INPUT  T$ 

877  DISP  9  DISP  "Whoro  do  you  want  lino": II  9  DISP  "prin 

tod?"  j 

878  OISP  " CGivo  y— coordinate  from  the"  9  DISP  "graph. ) "  j 

879  INPUT  T2  “  j 

880  MOVE  U1*U2*.  5,  T2  | 

881  LABEL  T$  | 

882  NEXT  II 

885  DISP  "Would  you  liko  another  copy?"  8  DISP  " CY/N) " 

890  INPUT  Q* 

900  IF  Q$-"Y"  THEN  685  j 

910  CLEAR  j 

920  OISP  "Are  you  finished?  CY/N)"  .  j 

930  INPUT  0$  / 

940  IF  Q$-"Y"  THEN  CLEAR  9  GOTO  960  j 
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I 


950 

960 

965 

970 

975 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1900 

1910 

1920 

1930 

1940 

1950 

1960 


GOTO  390 

OISP  "Would  you  lika  to  digitize" 
DISP  "another  set  of  curves?  (Y/N) " 
INPUT  Q$ 

IF  G$="Y"  THEN  120 
CHAIN  "SMILE:  0700" 

END 
REM  * 

REM  * 

REM  *  DIGITIZER  SUBROUTINE 
REM  *  A)  CHECK  STATUS: 

REM  *  B)  GET  POINTS:  OR. 

REM  *  C>  GET  SOFTKEYS. 

REM  * 

REM  * 

OUTPUT  706  :  "SG; OP" 

B-7  9  B$“"RS"  9  Bl-1 
P$«"BP24.  125.  5:  BP36.  125.  5" 

GOSUB  1200  9  REM  *  SOFTKEY? 

IF  S  THEN  1160 
B-2  9  B$-”OD"  9  Bl-3 
PS-"BP36.  125.  S" 

GOSUB  1200  9  REM  *  POINT? 

RETURN 
REM  * 

REM  * 

REM  *  STATUS  SUBROUTINE 
REM  *  A)  CHECK  SOFTKEY: 

REM  *  B)  CHECK  POINT. 

REM  * 

REM  * 

OUTPUT  706  :  "OS"  9  REM  *  READY? 
ENTER  706  :  S 
S-BITCS.  B) 

IF  NOT  S  THEN  GOTO  1370  9  REM  *  NO. 
OUTPUT  706  :  BS  9  REM  *  YES.  GET. 

IF  Bl-3  THEN  1350 
ENTER  706  :  TCI) 

GOTO  1360 

ENTER  706  ;  T  <15  .  T  C2>  .  T  C3) 

OUTPUT  706  : P$ 

RETURN 
REM  * 

REM  * 

REM  *  ACTIVE  SUBROUTINE 

REM  *  A)  ESTABLISH  ACTIVE  FILES: 

REM  *  S)  STORE  DATA. 

REM  * 

REM  * 


EXIT. 
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1970  GCSUB  1000  9  REM  *  DATA? 

1980  IF  NOT  S  THEN  GOTO  1970  9  REM  *  NO.  WAIT. 

1990  IF  B-2  THEN  GOTO  2030  9  REM  *  YES. 

2000  K-TC1) 

2010  ON  K  GOSUB  3100.3200.3300.3520.  3600.5900.  3800.3800. 

3800.  3800.  3800.  3800.  3800.  3800.  3800.  3800 

2020  GOTO  2070 

2030  GOSUB  4700 

2040  X  CPI.  P  CJ)  )  -T  Cl) 

2050  YCP1.  PCJ)  )  -TC2) 

2060  PCJ)-PCJ)+1 
2070  RETURN 
2100  REM  * 

21 1G  REM  * 

2120  REM  *  FILE  SUBROUTINE 

2130  REM  *  A)  CREATE  DATA  FILES. 

2140  REM  *  B)  OPENS  DATA  FILES. 

2150  REM  * 

2160  REM  * 

2170  CREATE  '‘X'-SVAL*  CQ)  .  150*0+10.  8 
2180  CREATE  " Y"&VAL$ CQ)  .  150*Q.  8 
2190  ASSIGN#  1  TO  "X"&VAL$CQ) 

2200  ASSIGN#  2  TO  "V&VALS  CQ) 

2210  RETURN 
2700  REM  * 

2710  REM  * 

2720  REM  *  MASS  STORAGE  SUBROUTINE 

2730  REM  *  A)  ESTABLISH  MEDIUM 

2740  REM  * 

2750  IF  R$*»"TAPE"  THEN  2780 
2755  IF  R$-,,0ISK01"  THEN  2770 
2760  MASS  STORAGE  IS  0700" 

2765  GOTO  2790 

2770  MASS  STORAGE  IS  ">0701 " 

2775  GOTO  2790 

2780  MASS  STORAGE  IS  ”« T" 

2790  RETURN 
2800  REM  * 

2810  REM  * 

2320  REM  *  STORAGE  SUBROUTINE 

2830  REM  *  A)  TRANSFER  FROM  ACTIVE 
2840  REM  *  FILE  TO  STORAGE  MEDIUM 

2850  REM  *  CT APE  OR  DISK) 

2860  REM  * 

2870  REM  * 

2880  FOR  1 1-1  TO  10 

2890  PRINT#  1.  II  :  PCI1) 

2900  NEXT  II  9  M-10 
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2901  FOR  14-1  TO  4 

2902  PRINT#  2.2*14-1  ,  U<I4) 

2903  PRINT#  2.2*14  j  VCI4) 

2904  NEXT  14 

2910  FOR  12-1  TO  PI 
2920  FOR  13-1  TO  PCI2)-1 
2930  M-M-l 

2940  PRINT#  l.M  s  XCI2.  13) 

2950  PRINT#  2.  M-2  j  Y«2.  13) 

2960  NEXT  13 
2970  NEXT  12 
2980  RETURN 
3000  REM  * 

3010  REM  * 

3020  REM  *  SOFTKEY  SUBROUTINE 

3030  REM  *  1)  STOP 

3040  REM  *  2)  STORE 

3050  REM  *  3)  LOAD 

3060  REM  *  4)  SCALE 

3070  REM  *  5)  PLOT 

3080  REM  *  6-16)  Not  used 

3090  REM  * 

3100  REM  * 

3110  REM  *  1)  STOP 

3120  REM  * 

3130  REM  * 

3140  CLEAR 

3150  DISP  "Are  you  finished  with  this  curve?"  3  DISP  " <Y 
/N)  '* 

3160  INPUT  AS 

3170  IF  AS— " Y“  THEN  GOTO  370 
3180  RETURN 
3190  REM  * 

3200  REM  * 

3210  REM  *  2)  STORE 

3220  REM  * 

3230  REM  * 

3240  CLEAR 

3250  OISP  "Where  do  you  want  the  data"  §  DISP  "stored?  < 
TAPE/DISKOO/OISKOl) " 

3260  INPUT  RS 
3270  GOSUB  2700 
3273  GOSUB  2100 
3280  GOSUB  2800 
3283  GOSUB  4200 
3286  GOSUB  5100 
3290  RETURN 
3300  REM  * 


3260  PRINT  "is  To-"» 

3262  YG-.  I*  CY  Cl.  P  Cl)  -1)  -V  C3)  )  +V  (3) 

3264  M0-CYC1.  I1  +  D-YC1.  I1))/CXC1.  Il-D-XCl,  ID) 

3265  IF  YO-YO,  115  <0  THEN  YC-Y(1.  ID 

3266  X0-CY0-YC1.  1 1)  ) /MO-OC  Cl.  ID 
3268  R1-CXQ-UCD)*N 

3270  GOSUB  6500 

3275  PRINT  VAL$CRl)j"  "j  U2$«  "  9  PRINT 

3277  PRINT  "Tha  value  of  instantaneous" 

3279  PRINT  V0$  Cl.  Z1  <15  3  i  “  at  To  is  "» 

3280  Rl-  CYO-V  C3)  )  *V1*V2  Cl)  /  CV  Cl)  -V  C3)  ) 

3291  GOSUB  6500 

3283  PRINT  VALSCRDj"  "j  V4$  Cl.  Z  Cl)  ]  ,  ”  9  PRINT 

3285  FOR  12-2  TO  4  STEP  2 
3287  13-12 

3289  IF  12-4  THEN  13-6 

3290  0-Vl*V2  CI2) / CV  Cl) -V  C3) ) * Cl+IP C CI2-1) /2) ) 

3292  MO—  CY  Cl  3,  I1  +  D-YCI3.  I1))/CXCI3.  I1  +  D-XCI3.  ID) 

3294  Y0-CX0-XC1.  1 1)  )  *MO-*-Y  CI3.  ID 

3295  PRINT  “Tha  value  of  instantaneous" 

3300  PRINT  V0SCZ1  CI2-D+1.  21  CI2)]»  “  at  To  is  "i 

3305  IF  12*2  THEN  3325 

3310  Rl- CYO-V C3))*0 

3315  GOSUB  6500 

3320  GOTO  3335 

3325  Rl -Y0*0/CK2#CUC2)-U  CD  >)  *N2 
3330  GOSUB  6500 

3335  PRINT  VAL$CRl)i"  "i  V4$  CZ  CI2-1)  +1.  Z  CI2)  ]  »  " 

3340  PRINT 
3345  NEXT  12 

3350  PRINT  "The  change  in  "i  V0$  CZ1  C3>  *1.  Z1  C4)  ]  ,  "  from  " 
3355  PRINT  "To  to  second— breakdown  is  " 

3360  R 1  -  CY  C6.  P  C6)  - 1)  - YO)  *0/  CK2*  CU  C2)  -U  C 1 )  )  )  *N2 
3365  GOSUB  6500 

3370  PRINT  VAL*CR1)j"  "i  V4$  CZ  C3)  +1.  Z  C4)  3  ,  "  9  PRINT 

3372  PRINT  "To  to  second-breakdown  is  " 

3374  R1-CXC1.  PCD-1)  -X0)  *N 

3376  GOSUB  6500 

3377  PRINT  VALSCRDi"  "j  U2*»  " 

3378  RETURN 

3380  CLEAR  9  GCLEAR  9  OISP  "Load  the  printer  and  press" 
3382  OISP  "CONTINUE. "  9  PAUSE 

3385  LORG  6  9  CSIZE  CH* CU C2) -U Cl ) ) /8831 *50+ 1 C) /19.  .  3 
3390  J-ia+IPCCIl  +  l)/4)*60 
3395  L-10+FP  Cl 1/2) *100 

3400  LOCATE  J.  H*  CU  C2)  -U  Cl)  )  /8831  *50-<-J.  L.  H*  CV  Cl) -V  C3)  )  /88 

31*50+L 

3405  FRAME 

3410  SCALE  0.  U1-U2.  0.  V1-V2CI1) 

3415  IF  G$-"Y"  THEN  3430 
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3015  R3-FNA CV2 C3) ) 

3020  N1-V2C3)/CR3*1000‘J5)  a  V2  C3)  -R3 
3025  K2-IP  C  <V  Cl)  -V  C3)  )  *2*  <U  (2?  -U  Cl)  )  /B) 

3030  IF  K2>0  THEN  K2-1/K2  9  GOTO  3040 
3035  K2-IPCB/CCVCD-VC3))  *2*(U(2)-U(1)))) 

3040  V2  C4) -K2*V2  Cl) *V2 <2) *V1*U1*U2 
3045  R4-FNACV2C4)) 

3050  N2-V2C4)/CR4*1000*J5)  9  V2  C4) -R4 
3055  J4C4)-J4C4)+J4C3)-2 
3060  G0SU8  3600 

3065  CLEAR  a  OISP  “What  Is  tha  horizontal  label?" 

3070  INPUT  US 

3075  OISP  a  OISP  "What  is  the  first  vortical" 

3080  DISP  "label?" 

3085  INPUT  VOS 

3090  OISP  a  DISP  "What  is  the  second  vertical" 

3095  OISP  "label?" 

3100  INPUT  VIS 

3105  21  CO)  “0  9  21  CD-LENCVOS) 

3110  V2S-" POWER"  a  V3S“" ENERGY" 

3115  21  C2)-LENCV1S)+21  Cl) 

3120  21  C3)-LENCV2S)+21  C2) 

3125  21  C4)-LENCV3$)>21  C3) 

3130  V0S-V0S&V1S&V2S&V3S 
3135  N-U1*U2/CUC2)-UC1)) 

3140  FOR  1 1-1  TO  4 

3145  0— V 1  *V2  CID/CVCl)  —V  C3)  )  *  Cl  +  IP  C  Cl  1-1)  /2)  ) 

3150  PRINT  "Tha  value  for  instantaneous" 

3155  PRINT  V0SC21  CI1-D+1.  21  CIl)3j  "  at  the  final  point 

s  " 

3160  IF  I 1 >2  THEN  3180 

3165  Rl  —  CY  <11.  PCD— 1)  -V  C3)  )  *0 

3170  GOSUB  6500 

3175  GOTO  3210 

3180  IF  II >3  THEN  3200 

3185  R1-YC5.  PC5)-1)*0/K1-*N1 

3190  GOSUB  6500 

3195  GOTO  3210 

3200  Rl— Y  C6.  P  C6)  -1)  *0/  CK2*  CU  C2)  -U  Cl)  )  )  *N2 
3205  GOSUB  6500 

3210  PRINT  VALSCR1),"  "j  V4S  C2  Cl  l-l)  *1.  2  Cl  1)  3 1  "  at  ", 
3215  R1-CXC1.  PC1)-1)-UC1))*N 
3220  GOSUB  6500 

3225  PRINT  VALSCRD,"  ",  U2S,  ".  "  9  PRINT 

3230  NEXT  II 

3235  FOR  11-2  TO  PC1)-1 

3240  IF  CYC1.  ID -VC3))/CYC1.  P  Cl) -1) -V  C3)  )  >.  1  THEN  3247 
3245  NEXT  II 
3247  Il-Il-l 

3250  PRINT  a  PRINT  "The  time  at  10X  " 
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2730 

REM  *  1) 

MAKES  SEPARATE 

PLOTS 

2740 

REM  * 

OF  2  OBJECT  CURVES. 

2750 

REM  * 

PRODUCT  CURVE. 

AND 

2760 

REM  * 

TIME  INTEGRAL 

2770 

REM  * 

CURVE 

2780 

REM  * 

2790 

REM  * 

2800 

11-1 

2805 

GOTO  3380 

2810 

11-2 

2815 

GOTO  3380 

2820 

11-3 

2825 

GOTO  3380 

2830 

11-4 

2835 

GOTO  3380 

284C  CLEAR  9  GCLEAR 
2842  8EEP  200.  50 

2845  DISP  “Do  you  want  a  standard  8  x  10" 

2850  OISP  “grid?  (Y/N) “ 

2855  INPUT  Gl$ 

2860  IF  G1$-"Y"  THEN  G9-"Y"  9  Ul-10  9  Vl-3  9  GOTO  2900 
2865  CLEAR  9  DISP  "Do  you  want  a  grid  on  the  plot?" 
2870  DISP  " (Y/N) " 

2875  INPUT  CS 

2880  DISP  9  DISP  "How  many  horizontal  division®?'* 

2885  INPUT  U1 

2890  DISP  9  DISP  "How  many  vortical  divisions?" 

2895  INPUT  VI 

2900  CLEAR  9  DISP  "How  much  is  each  horizontal" 

2905  DISP  "division?  (magnitude,  unit) " 

2910  DISP  "  CS  mV  is  5.  mV.  3  " 

2915  INPUT  U2.  U2S 

2920  DISP  9  OISP  "How  much  is  each  division  on” 

2925  DISP  "tha  first  vertical  axis?" 

2930  DISP  “  (magni tudo,  uni t) " 

2935  DISP  " C5  mV  is  5.  mV.  3 " 

2940  INPUT  V2(1).V4$ 

2945  DISP  9  DISP  "How  much  is  each  division  on" 

2950  OISP  "the  second  vertical  axis?" 

2955  OISP  "  (magnitude,  uni t)  ” 

2960  DISP  “  C5  mV  is  5.  mV.  3  " 

2965  INPUT  V2(23.V5* 

2970  Z(0J-0  9  Z(1)-LEN(V4*> 

2975  Z (2)  -LEN (VS*) *Z  (l) 

2980  V4*-V4*&V5* 

2985  Ml “0  9  Ll-0  9  J-3 
2990  MS-"  "  9  L*-"  "  9  J*-"  " 

2995  K1-IP((V(1)-V  (3)  )  *2/ A) 

3000  IF  K1>0  THEN  Kl-l/Kl  9  GOTO  3010 
3005  K1-IP(A/(V(1)-V(3>)  *23 
3010  V2(3)-K1*V2(13  wV2(23*Vl 


-78- 


1410  IF  YCA1.  Il)>-A  THEN  A-YCAl.  ID 

1420  NEXT  II 

1440  FOR  11-1  TO  PC15-1 

1450  X  <81.  1 1  >  -X  Cl.  ID 

1470  IF  1 1  —  1  THEN  YCBl.  I15-YCA1.  I1)*(X<1,  Il)-U (3))  3  GOT 
0  1490 

1480  YCBl.  I15-YCA1.  I1)*CX<1.  ID-XC1.  1 1-D  )  +Y  <B1.  Il-D 

1485  IF  XC1.  ID-XC1.  Il-1)>-C  THEN  OXCl.  I15-XC1.  Il-D 

1490  IF  YCBl.  Il)*-B  THEN  B-Y  CB1.  ID 

1500  NEXT  II 

1510  RETURN 

1600  REM  * 

1610  REM  * 

1620  REM  *  READ  SUBROUTINE 
1630  REM  *  D  READS  ARRAY 
1640  REM  *  CTAPE  OR  DISK) 

1650  REM  * 

1660  REM  * 

1670  FOR  12-1  TO  PI 
1680  FOR  13-1  TO  P <125-1 
1690  K-K-l 

1700  READ#  l.K  i  XCI2.  13) 

1710  READ#  2.  K-2  i  YCI2.  12) 

1720  NEXT  13 
1730  NEXT  12 
1740  RETURN 
1800  REM  * 

1810  REM  * 

1820  REM  *  CLOSE  SUBROUTINE 

1830  REM  *  1)  CLOSE  FILES 

1840  REM  * 

1850  REM  * 

1860  ASSIGN#  1  TO  * 

1870  ASSIGN#  2  TO  * 

1880  RETURN 
1900  REM  * 

1910  REM  * 

1920  REM  *  MASS  STORAGE  SUBROUTINE 
1930  REM  *  15  ESTABLISH  MEDIUM 

1940  REM  * 

1950  IF  R$-"TAPE"  THEN  1980 
1955  IF  R$-"DISK01 "  THEN  1970 
I960  MASS  STORAGE  IS  0700" 

1965  GOTO  1990 

1970  MASS  STORAGE  IS  ".0701” 

1975  GOTO  1990 

1980  MASS  STORAGE  IS  T” 

1990  RETURN 
2700  REM  * 

2710  REM  * 

2720  REM  *  PLOT  SUBROUTINE 


1 130  X  <4,  ID  *X  Cl.  ID 

1135  NEXT  II 

1140  FOR  1 1-1  TO  P<25-1 

1145  X  <4.  1 1*P  <15 -15 -X  <2.  ID 

1150  NEXT  II 

1155  FOR  1 1-1  TO  P<45-1 

1160  FOR  12-11  TO  P<45-1 

1165  IF  X<4.  I15>X<4.  125  THEN  R5-X  C4,  ID  9  X  C4,  1 15 -X  C4.  12 
5  0  X  (4.  125 -R5 
1170  NEXT  12 
NEXT  II 

FOR  11-2  TO  P  <15 -1 
FOR  12-14  TO  P  <45 -2 
14-12 

IF  X  <4.  125  >X  <1.  1 15  THEN  1210 
Y<4.  12-15  — M  <1.  11  —  15  *X  <4,  125  *M  <2.  11-15 
NEXT  12 
NEXT  II 

P  <15  -P  <45  —2  Q  14-2 
FOR  1 1  —  1  TO  P  <15 -1 
Y<1.  1 1 5  —  Y  <4.  ID 
X<1.  1 15  —X  <4.  11*15 
NEXT  II 

FOR  11-2  TO  P  <25 -1 
FOR  12-14  TO  P  <45 -2 
14-12 

IF  X  <4.  125  >X  <2.  1 15  THEN  1255 

Y  <4.  12-15  -M  <3.  1 1-15  *X  <4.  125  *M  <4.  11-15 
NEXT  12 
NEXT  II 

P  <25  — P  <45  —2  9  P  <35  -P  <25 
FOR  1 1-1  TO  P  <25 -1 

Y  <2.  1 15  — Y  <4.  ID 
X  <2.  1 15  —X  <4,  11*15 
NEXT  II 

FOR  1 1-1  TO  P  <35 -1 
X  <3.  1 15  —  Y  <C2.  ID 

Y  <3.  1 15  — Y  <C1,  ID 
NEXT  II 


1175 

1180 

1165 

1190 

1195 

1200 

1205 

1210 

1215 

1217 

1219 

1221 

1223 

1225 

1230 

1235 

1240 

1245 

1250 

1255 

1260 

1262 

1264 

1266 

1268 

1270 

1275 

1280 

1285 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 


SUBROUTINE 

EVALUATE  PRODUCT  AND 
TIME  INTEGRAL  FOR  CURVES 


RETURN 
REM  * 

REM  * 

REM  *  MATH 
REM  *  15 

REM  * 

REM  * 

REM  * 

a— o  a  c-o 

FOR  11-1  TO  P  <1 5 -1 
X  <A1 ,  1 15  —X  Cl.  ID 

Y<A1.  ID  —  <Y  <1,  1 1 5  -V  <35  5  *  <Y  <2.  1 15 -V  <35  5 
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445  GOSUB  5300 

450  Al-5  9  Bl-6 

490  GOSUB  1000 

510  p(5)«p(d  a  p<65-p<15 

520  GOSUB  1300 

525  GOSUB  2B40 

530  ON  KEY#  1.  "PLOT  V"  GOTO  2B00 

535  ON  KEY#  2.  "PLOT  I"  GOTO  2B10 

540  ON  KEY#  3.  "PLOT  P"  GOTO  2020 

545  ON  KEY#  4.  "PLOT  U"  GOTO  2B30 

550  ON  KEY#  5.  "V-I"  GOTO  5600 
555  ON  KEY#  6. "SCALE"  GOTO  4400 
560  ON  KEY#  7. "STOP"  GOTO  600 
570  CLEAR  a  KEY  LABEL 
575  BEEP  200. 50 

500  DISP  "Press  the  desired  key.  '* 

590  GOTO  590 

600  DISP  “Are  you  finished  with  "  a  OISP  "this  operation 
7  CY/N5  " 

610  INPUT  AS 

620  IF  AS#" Y"  THEN  530 

710  DISP  "Do  you  want  to  work  on" 

720  DISP  "another  set  of  curves?" 

730  DISP  " CY/N5 " 

740  INPUT  AS 

750  IF  AS-"Y"  THEN  10 

755  CLEAR  8  MASS  STORAGE  IS  ".0700" 

760  CHAIN  "SMILE" 

770  END 
1000  REM  * 

1010  REM  * 

1020  REM  *  INTERPOLATION  SUBROUTINE 
1030  REM  *  1)  FIND  SLOPES  AND  INTERCEPTS 

1040  REM  *  2)  EVALUATE  SLOW  GRAPH 

1050  REM  *  FOR  RAPID  GRAPH  POINTS 

1060  REM  • 

1070  REM  * 

1080  FOR  1 1-1  TO  P  <15 -2 

1082  IF  XC1.  I1-15-XC1.  115-0  THEN  M  Cl.  1 15 -M  <1.  1 1-15  8  GOT 
0  1090 

1085  MCI.  I15-CYC1.  I1  +  15-YC1.  I155/CXC1.  I1+15-XC1.  1155 

1090  M  <2.  115 -Y  Cl.  115  -M  Cl.  I15*XC1.  115 

1095  NEXT  II 

1100  FOR  1 1-1  TO  PC25-2 

1102  IF  X  <2.  Il-*15— XC2.  115—0  THEN  M  <3,  1 15  -M  <3.  1 1-15  9  GOT 
0  1110 

1105  MC3.  I15-CYC2,  I1  +  15-YC2.  I15  5/CXC2.  I1-15-XC2.  115  5 
1110  MC4.  I15-YC2.  I15-MC3.  I15*XC2.  115 
1115  NEXT  II 

1120  14-2  9  PC45-PC15+PC25-1 
1125  FOR  1 1-1  TO  PC15-1 
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10  REM  * 

20  REM  * 

30  REM  *  MATH  PACKAGE 
40  REM  *  MAIN  PROGRAM 
50  REM  * 

60  REM  *Copyrightt  11/16/93 
70  REM  *  gandalf  software,  inc. 

SO  REM  *  Chuck  Graves,  wizard 
90  REM  * 

100  REM  * 

105  SHORT  X  <6.  1 50)  .  Y  (6.  1 50)  .  M  (4.  1 50) 

110  INTEGER  P  (10)  .  U  (4)  .  V  (4)  ,  Z  (4)  .  Z1  C4)  .  J4  (4) 

115  DIM  VO$C303 
120  COSUB  4500 
125  BEEP  200.50 

130  DISP  "What  is  the  name  of  the" 

140  OISP  "x-file?" 

150  INPUT  XS 

160  IF  LEN(XS)»10  THEN  DISP  "Name  is  too  large." 
130 

180  OISP  9  OISP  "What  is  the  name  of  the" 

190  OISP  "y-file?" 

200  INPUT  Y$ 

210  IF  LEN  <YS) >10  THEN  DISP  "Name  is  too  large." 
180 

220  DISP  3  OISP  "Who re  ore  "»  XS,  "  and  "i 
230  OISP  Y$,  "  stored?  (TAPE/DISKOO/DISKOl) “ 

240  INPUT  RS 

250  GOSUB  1900 

260  ASSIGN#  1  TO  X* 

270  ASSIGN#  2  TO  Y$ 

280  FOR  1 1-1  TO  10 
290  READ#  1.  II  i  PCI1) 

300  IF  P(I1)»1  THEN  Ql-Il  8  Pl-Il 

310  NEXT  II  a  K-10 

320  FOR  12-1  TO  4 

330  READ#  2.2*12-1  i  U<12? 

340  READ#  2.  2*12  i  V<12) 

350  NEXT  12 
360  GOSUB  1600 
370  GOSUB  1800 
380  CLEAR 
385  BEEP  200.  50 

390  DISP  "Please  list  the  two  curves  you” 

400  OISP  "wish  multiplied.  (Separate  the" 

410  OISP  "curves  with  a  comma  and  in  the"  8  DIS1 
VOLTAGE.  CURRENT.  > " 

420  INPUT  A.  8 
425  CLEAR 

440  OISP  "WORKING" 


8  GOTO 


8  GOTO 


"or  den 
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6135  DISP  9  OISP  "Enter  forward  base  currant.  " 

6140  INPUT  0$a  PRINT  0* 

6145  PRINT  a  PRINT  "Reverse  baaa  current!  " 

6150  PRINT  a  PRINT  "  Nominal  -  “» 

6155  OISP  8  OISP  "Entar  nominal  reverse  base  currant 
6160  INPUT  0*0  PRINT  0$ 

6170  OISP  a  OISP  "Enter  ravarsa  base  currant  at  SB" 
6180  INPUT  0*8  PRINT  "  Actual  at  SB  -  "j  D* 

6190  PRINT  a  PRINT  "Comments!  " 

6195  DISP  3  OISP  "Enter  any  additional  comments.  “ 
6200  INPUT  0*3  PRINT  ”  ",  0$ 

6205  PRINT  a  PRINT  3  PRINT 
6210  RETURN 


5760 

5770 

S780 

5790 

5800 

5810 

5820 

5830 

5840 

5841 

5842 

5843 
58S0 
5860 
5870 
5880 
5885 
5890 
5900 
5910 
5920 
5930 
5940 
5950 
5960 
5965 
5970 
5975 
5980 
5990 
5995 
6000 
6010 
6015 
6020 
6025 
6030 
6040 
6050 
6055 
6060 
6070 
6075 
6080 
6090 
6095 
6100 
6105 
6110 
6120 
6130 


PRINT  9  PRINT  "3  LOAD' 


PRINT 

PRINT 


PRINT  a  PRINT 
PRINT  "  Cha 
PRINT  ”  of 


PRINT 

PRINT 


PRINT  a  PRINT 


Load  old  curves" 
from  tape  or  disk. " 


SCALE" 


Change  scale  size' 
of  plot.  " 


PRINT  3  PRINT  "5  PLOT* 


Make  intermedi ate” 
plot  of  curves. “ 


DEVICE  DATA" 


PRINT  "  Make  printout  of" 

PRINT  "  device  data.  " 

PRINT  8  PRINT  "7-16  Not  used.  " 

FOR  1 1-1  TO  5 
PRINT  a  PRINT 
NEXT  II 
CLEAR 
RETURN 
REM  * 

REM  * 

REM  *  DEVICE  DATA 
REM  -  SUBROUTINE 
REM  * 

REM  * 

CLEAR 

PRINTER  IS  2 

PRINT  a  PRINT  "Manufacturer  -  "i 
DISP  "Enter  manufacturer. " 

INPUT  0*8  PRINT  0$ 

PRINT  8  PRINT  "Oevice  type  -  ". 

DISP  a  DISP  "Enter  device  type.  " 

INPUT  0*8  PRINT  0* 

PRINT  a  PRINT  "Mask  type  or  other  applicable 
DISP  8  OISP  "Enter  mask  type  or  other" 

PRINT  "data  -  ", 

OISP  "applicable  data." 

INPUT  0*8  PRINT  D* 

CLEAR 

PRINT  a  PRINT  "Device  number  -  "» 

OISP  a  DISP  "Enter  device  number.  " 

INPUT  D$8  PRINT  0* 

PRINT  a  PRINT  "Temperature  -  "» 

DISP  8  OISP  "Enter  temperature  CRT  or" 

DISP  "degrees  C5  .  " 

INPUT  0* 

PRINT  0*s "  C" 

PRINT  8  PRINT  "SB  type  -  "i 

DISP  a  OISP  "Enter  SB  type  CN. A. B  or  C) . " 

INPUT  D*a  PRINT  D* 

CLEAR 

PRINT  8  PRINT  "Forward  base  current  -  "i 
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5200  DISP  "What  name  do  you  want  for  the"  9  OISP  "y-coor 
dinate  file?" 

5210  INPUT  Y5$ 

5220  IF  LEN  <Y5$)  >10  THEN  DISP  "Noma  is  too  large."  9  GOT 
0  5200 

5230  RENAME  "X“&VALSCQ)  TO  X5* 

5240  RENAME  "Y"aVAL$CQ)  TO  Y5$ 

5250  RETURN 
5300  REM  * 

5310  REM  * 

5320  REM  *  INITIALIZATION  SUBROUTINE 
S330  REM  * 

5340  REM  * 

S34S  OISP  9  OISP  "INITIALIZING" 

5350  FOR  I 1-1  TO  Q 

5360  FOR  12-1  TO  150 

5370  Xdl.  12) -0 

5380  Y  <1 1.  12) -0 

5390  NEXT  12 

5400  NEXT  II 

5410  FOR  I 1-1  TO  4 

5420  UCID-0 

S43o  van -a 

5440  NEXT  II 

5450  FOR  I 1-1  TO  10 

S460  PCID-l 

5470  NEXT  11 

5480  Pl-1  9  Nl-0  9  H-l  9  H$-"A" 

5490  G1-. 0213  9  G2-. 022 
5500  PLOTTER  IS  1 
5505  CLEAR 
5510  RETURN 
5600  REM  * 

5610  REM  * 

5620  REM  *  VISUAL  AID  SUBROUTINE 

5630  REM  *  A)  PRINT  SOFTKEY  ASSIGNMENTS. 

5640  REM  * 

5650  REM  * 

5660  PRINTER  IS  2 

5670  PRINT  "Softkey  Assignments" 

5680  PRINT  9  PRINT 
5690  PRINT  "1  STOP" 

5700  PRINT  "  End  of  curve." 

5710  PRINT  "  Depress  to  start" 

5720  PRINT  "  next  curve.  " 

5730  PRINT  9  PRINT  "2  STORE” 

5740  PRINT  "  Transfer  curves" 

5750  PRINT  ”  to  tape  or  disk.  " 
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4310  REM  * 

4320  REM  *  FUNCTION  FNX 
4330  REM  * 

4340  REM  * 

4350  REM  * 

4360  DEF  FNX<A5 
4370  B-H*A 
4380  FNX-IP  <B-<-.  55 
4390  FN  ENO 
4500  REM  * 

4510  REM  * 

4520  REM  *  CORRECTION  SUBROUTINE 

4530  REM  *  A)  CALCULATE  CORNERS.  <UL.  UR.  LL.  LR5 

4540  REM  *  B5  CALCULATE  CORRECTION  FACTORS. 

4550  REM  * 

4560  REM  * 

4570  V 1  -SQR  <  <U  <35  — U  <45  5  ‘2*  <V  <35 -V  <45  5  *2) 

4580  US-ATN  <  <V  <45  -V  <35  5  /  <U  <45  -U  <35  5  5 
4590  U  <45  -U  <35  ♦  I P  CV  1  -*•.  5) 

4600  V<45-V<35 

4610  Vl-SQR  <  <U  <15  -U  <35  5  *2+  <V  <15  -V  <35  5  *2) 

4620  U(l)-UQ) 

4625  V  <  1 5  —V  <35  -*■  I P  <V  1  5) 

4630  U<25-U<45 
4640  V  <25  “V  <15 
4680  RETURN 
4700  REM  * 

4710  REM  * 

4720  REM  *  SKEW  SUBROUTINE 
4730  REM  *  A5  ADJUST  ROTATION 
4740  REM  * 

4750  REM  * 

4760  Vl-SQR  <  <U  <35  -T  <15  5  ‘2+  <V  <35  -T  <25  5  ‘25 

4765  IF  T <15 -U <35-0  THEN  V5-PI/2  9  GOTO  4780 

4770  V5— ATN  <  <T  <25  -V  <35  5  /  <T  <15  -U  <35  5  5 

4780  IF  V5<0  AND  T <25 -V <35 >0  THEN  V5-V5+PI 

4790  T  <  1 5  -U  <35  -*•  I P  <V  1  *COS  <V5-U55  * .  55 

4800  T  <25  -V  <35  +  IP  <V1*SIN  <V5-U55  55 

4810  RETURN 

5100  REM  * 

5110  REM  * 

5120  REM  *  RENAME  SUBROUTINE 

5130  REM  *  A5  REGISTERS  ARE  USER  NAMED. 

5140  REM  * 

5150  REM  * 

5160  CLEAR 

5170  DISP  "What  name  do  you  want  for  tha"  9  DISP  "x-coor 
dinata  file?" 

5180  INPUT  X5* 

5190  IF  LEN<X5*5>10  THEN  DISP  "Noma  is  too  large."  9  GOT 
0  5170 
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3730  ORAW  FNX  <X  <  1 1 .  1 25  5  .  FNX (Y  ( 1 1 .  12)) 

3740  NEXT  12 
3750  NEXT  II 

3752  MOVE  0.  V  <35 

3753  ORAW  12032.  V  (3) 

3755  WAIT  4000 
3760  RETURN 

3770  REM  * 

3780  REM  * 

3790  REM  *  6-16)  Not  used 

3800  REM  * 

3810  REM  * 

3820  CLEAR 

3830  DISP  9  OISP  "This  softkey  performs  no  function. 
3840  OISP  9  OISP  "Please  try  another. 

3850  Nl-Nl+1 

3860  IF  N1 <4  THEN  GOTO  1900 
3865  PLOTTER  IS  1 
3870  CSIZE  21  9  GCLEAR 
3880  SCALE  0.  10.  0.  10 
3890  MOVE  1. 5 
3900  LABEL  "HANDS" 

3910  MOVE  2. 0 
3920  LABEL  "OFF" 

3930  WAIT  2000  9  GCLEAR 
3940  RETURN 
4000  REM  * 

4010  REM  * 

4020  REM  *  READ  SUBROUTINE 
4030  REM  *  1)  READS  ARRAY 

4040  REM  *  <TAPE  OR  DISK) 

4050  REM  * 

4060  REM  * 

4070  FOR  12-1  TO  PI 
4080  FOR  13-1  TO  P<12)-1 
4090  M-M-l 

4100  READ#  l.M  i  XCI2.  13) 

4110  READ#  2.  M-2  ;  YCI2.  13) 

4120  NEXT  13 
4130  NEXT  12 
4140  RETURN 
4200  REM  * 

4210  REM  * 

4220  REM  *  CLOSE  SUBROUTINE 

4230  REM  *  1)  CLOSE  FILES 

4240  REM  * 

4250  REM  * 

4260  ASSIGN#  1  TO  * 

4270  ASSIGN#  2  TO  * 

4280  RETURN 
4300  REM  * 
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3)  LOAD 


3320  REM  * 

3330  REM  * 

3340  REM  * 

3350  CLEAR 

3360  DISP  "What  is  the  namo  af  the  x-file?" 

3370  INPUT  X$ 

3375  IF  LENCX$)>10  THEN  DISP  “Name  is  too  larqe. "  9  GOTO 
3360 

3380  OISP  i  DISP  "Whet  is  ths  nama  of  the  y-file?" 

3390  INPUT  Y* 

3395  IF  LENCY$)>10  THEN  DISP  "Namo  is  too  larqe. "  9  GOTO 
3380 

3400  DISP  9  DISP  "Whore  are  X$;  "  and  "i  9  DISP  Y$,  "  sto 
rod?  CTAPE/OISKOO/DISKOU  " 

3410  INPUT  RS 

3420  GOSUB  2700 

3430  ASSIGN#  1  TO  X$ 

3440  ASSIGN#  2  TO  Y$ 

3450  FOR  1 1-1  TO  10 

3460  READ#  1.  II  ,  PCI  13 

3470  IF  PCI1)»1  THEN  Q-Il  9  Pl-Il 

3480  NEXT  II  9  M-10 

3482  FOR  12-1  TO  4 

3484  READ#  2.2*12-1  ,  UCI25 

3486  READ#  2.2*12  i  VCI25 

3488  NEXT  12 

3500  GOSUB  4000 

3505  GOSUB  4200 

3510  GOTO  390 

3520  REM  * 

3530  REM  * 

3540  REM  *  4)  SCALE 

3550  REM  * 

3560  REM  * 

3570  CLEAR 

3580  OISP  "Enter  the  scale  factor  for"  9  DISP  "plotting. 

Close  than  4.  0)  " 

3590  INPUT  H 
3600  REM  * 

3610  REM  * 

3620  REM  *  5)  PLOT 

3630  REM  * 

3640  REM  * 

3650  CLEAR 

3660  OISP  "Do  you  want  a  plot  or  display"  9  DTSP  "of  the 
data  so  far?"  9  OISP  " CD  I SP/PLQT) " 

3670  INPUT  H$9  CLEAR  9  GCLEAR 

3690  SCALE  0.12032.0.9012 

3700  FOR  1 1-1  TO  PI 

3710  MOVE  FNXCXCI1.  1)).  FNXCYCI1.  1)) 

3720  FOR  12-2  TO  PCI1)-1 
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3420  AXES  U2.  V2CI1).  0.  0.  1.1.3 
3425  GOTO  3435 

3430  GRID  U2.  V2CI1).  0.  0.  1.  1.  1 

3435  FOR  12-1  TO  U1 

3440  MOVE  I2*U2.  -CV2CI1)  /10) 

3445  LABEL  VAL$CI2*U2) 

3450  NEXT  12 

3455  LORG  8 

3460  FOR  13-1  TO  VI 

3465  MOVE  -  CU2/8) .  I3*V2 Cl 1) 

3470  LABEL  VAL$ CIP <1 3*V2 Cl 1 ) ■*. 5) ) 

3475  NEXT  13 

3480  LORG  6  9  CSI2E  CH*  CU  <25  -U  Cl)  )  /8831*100-«-10) /16.  .  5 
3485  MOVE  . 5-Ul*U2. - CV1*V2 Cl  1) /10) 

3490  LABEL  U$,  "  C",  U29i  ")  " 

3495  DEG  9  LDIR  90  9  LORG  4 
3500  MOVE  -CU1*U2/10).  .  5*VI*V2CI1) 

3505  LABEL  V0$ CZ1 Cl  1-1) +1. 21  (1 15 3 j  "  C", V4* CZ <1 1 -1) +  1. Z Cl 
1)3 i ") M 

3510  RAD  9  LDIR  0 

3515  0-Vl*V2CIl)/CVCl)-VC3))  *  Cl-<-IP<CIl-l)/2)) 

3520  IF  11*2  THEN  3550 

3525  MOVE  <X  Cl  1.  1)  -U  Cl)  )  *N.  CY  Cl  1.  1) -V  C3>  )  *0 
3530  FOR  12-2  TO  PCI1)-1 

3535  DRAW  CXCI1.  I2)-UC1))*N.  CYCI1.  I2)-VC3))*0 

3540  NEXT  12 

3545  GOTO  3600 

3550  IF  I 1*3  THEN  3580 

3555  MOVE  CX  <5.  1)  -U  <1)  )  *N.  Y  <5.  1 )  *0/Kl *N1 
3560  FOR  12-2  TO  PCS) -l 

3565  DRAW  CX  <5.  1 2)  -U  C 1 )  )  *N.  Y  C5.  1 2)  *0/K  1  *N  1 
3570  NEXT  12 
3575  GOTO  3600 

3580  MOVE  CXC6.  1)-UC1))*N.  Y  C6.  1)  *0/ CK2*  CU  C2) -U  Cl)  )  )  *N2 
3585  FOR  12-2  TO  PC6)-i 

3590  DRAW  CX  C6.  12) -U  Cl)  )  *N.  Y  C6.  12)  *0/ CK2*  CU  C2) -U  Cl)  )  )  *N2 
3595  NEXT  12 

3600  CLEAR  9  DISP  "How  many  linos  doss  the" 

3605  DISP  "main  titls  havo?" 

3610  BEEP  200. 50 
3615  INPUT  T1 
3620  LORG  4 
3625  FOR  12-1  TO  T1 

3630  DISP  9  DISP  "What  is  lino  ",  VAL$ CI2) ,  "?" 

3635  INPUT  T* 

364C  DISP  9  DISP  "Where  do  you  want  lino  ",VALSCI2) 

3645  DISP  "printed?" 

3650  DISP  "CGlve  y-coordinate  from  the" 

3655  DISP  "graph.)" 
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3660  INPUT  T2 

3665  MOVE  . 5*U1*U2.  T2 

3670  LABEL  T$ 

3675  NEXT  12 
3680  GOTO  530 
3800  REM  •* 

3810  REM  * 

3820  REM  *  STRING  SUBROUTINE 
3830  REM  *  1)  CALCULATE  UNITS  FOR 

3840  REM  *  PRODUCT  AND  INTEGRAL 

3850  REM  *  CURVES 

3860  REM  * 

3870  REM  * 

3875  IF  L£NCJ$)<2  THEN  3890 
3885  L$~L$&J$  CJ4  <2> .  J4  C3) -13 
3890  IF  ZC1)<-1  THEN  3910 
3900  L$-L«V4$C1.2a)-l] 

3910  IF  ZC2)-ZC1)<-1  THEN  3930 
3920  L$-L*&V4$  CZ  Cl)  -"l.  Z  (2)  -13 
3930  IF  LEN<L$X-2  THEN  4060 
3940  FOR  1 1-1  TO  LEN <L$) 

3943  IF  L$  Cl 1.  113— "M"  THEN  Ll-Ll+6 
3946  IF  L*CI1.  113— "W  THEN  Ll-Ll+3 
3950  IF  L*  C 1 1 .  1 1 3  - " m ”  THEN  Ll-Ll-3 
3955  IF  LSCI1.  11.1 -"u"  THEN  Ll-Ll-6 
3960  IF  L$  Cl  1.  113  — “n"  THEN  Ll-Ll-9 
3970  IF  L*CI1.  1 1 3  —  •* p “  THEN  Ll-Ll-12 

3980  NEXT  II 

3981  IF  Ll-12  THEN  L*-"  T" 

3982  IF  LI -9  THEN  LS-"  G" 

3983  IF  LI -6  THEN  LS-"  M" 

3984  IF  LI -3  THEN  L*-M 

3985  IF  LI -0  THEN  LS-"  ** 

3986  IF  LI— 3  THEN  L$-“  m" 

3990  IF  LI— 6  THEN  L$-"  u” 

4000  IF  LI— 9  THEN  L$-"  n “ 

4010  IF  LI  — 12  THEN  L*-‘*  p" 

4020  IF  LI  — 15  THEN  L*-"  f" 

4030  IF  LI  — 18  THEN  L$-"  c” 

4040  IF  LI— 21  THEN  L*-"  y" 

4050  IF  LI— 24  THEN  L$-"  z" 

4060  L$-L«&”WM 

4065  IF  LENCJSX2  THEN  4075 
4070  MS-MS&J$ CJ4 C3) . J4 C4> -13 
4075  IF  ZC1X-1  THEN  4085 
4080  M*-M*&V4$C1.  ZCl) -13 
4085  IF  Z<2)-ZC1X-1  THEN  4092 
4090  M*-M$&V4*  CZ  <1X1.  Z  C2) -13 
4092  IF  LENOJ2SX-1  THEN  4100 
4095  MS-M5&U2* Cl.  LEN CU2*) -13 
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4100  IF  LENCMSX-2  THEN  4310 
4110  FOR  1 1-1  TO  LEN  CMS) 

4120  IF  MSCIl. I13--M"  THEN  Ml-Ml+6 
4130  IF  MSCIl.  I13-"k"  THEN  Ml-Ml+3 
4140  IF  MSCIl.  in-"m"  THEN  Ml-Ml-3 
4145  IF  MSCIl. I13-"u"  THEN  Ml-Ml-6 
4150  IF  MS  Cl 1.  1 13  — "n“  THEN  Ml-Ml-9 
4160  IF  MS  Cl 1.  I 13 -"p"  THEN  Ml-Ml-12 
4170  NEXT  II 

4180  IF  Ml-12  THEN  MS-"  T" 

4190  IF  Ml “9  THEN  MS-"  G" 

4200  IF  Ml -6  THEN  MS-"  M" 

4210  IF  Ml “3  THEN  MS-"  k" 

4220  IF  Ml-0  THEN  MS-"  " 

4230  IF  Ml— 3  THEN  MS-"  m" 

4240  IF  Ml— 6  THEN  MS-"  u" 

4250  IF  Ml— 9  THEN  MS-"  n" 

4260  IF  Ml  — 12  THEN  MS-"  p" 

4270  IF  Ml  — 15  THEN  MS-"  f" 

4280  IF  Ml  — 18  THEN  MS-“  a" 

4290  IF  Ml— 21  THEN  MS-"  y" 

4300  IF  Ml— 24  THEN  MS-"  z" 

4310  MS— MS&" J" 

4320  MS-MS  C2.  LEN  CMS)  3 
4330  LS-LS  C2.  LEN  CLS)  3 
4340  Z  C3)  -LEN  CLS)  +Z  C2) 

4350  Z  C4) -LEN  CMS) +Z  C3) 

4360  V4S-V4S&LS&MS 
4370  RETURN 
4400  REM  * 

4410  REM  * 

4420  REM  *  SCALE  SUBROUTINE 
4430  REM  * 

4440  REM  * 

4450  CLEAR 

4460  DISP  "Entar  tha  seals  factor  for" 

4470  DISP  "plotting.  Class  than  4.0)" 

4480  INPUT  H 

4485  H-2-»H 

4490  GOTO  530 

4500  REM  * 

4510  REM  * 

4520  REM  *  INITIALIZATION  SUBROUTINE 
4530  REM  * 

4540  REM  * 

4550  CLEAR  9  GCLEAR 

4560  DISP  8  DISP  "INITIALIZING" 

4565  FOR  1 1-1  TO  2 
4570  FOR  12-1  TO  150 
4575  XCI1.  12) -0 
4580  YCI1.  12) -0 
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4585  NEXT  12 
4590  NEXT  II 
4630  FOR  1 1-1  TO  4 
4635  Zai)-0 
4640  UCID-0 

4650  van -a 

4655  J4<I15-2 
4660  NEXT  II 
4670  FOR  1 1-1  TO  2 
4680  pan-i 
4690  NEXT  11 
4700  Pl-1  9  H— 1 
4710  PLOTTER  IS  705 
4725  FOR  11-3  TO  6 
4730  FOR  12-1  TO  150 
4735  MCI.  12)  -0 
4740  Xai.  125-0 
4745  M  (2,  125-0 
4750  Y  Cl  1.  125-0 
4755  NEXT  12 
4760  NEXT  II 
4765  FOR  11-3  TO  10 

4770  pan-i 

4775  NEXT  II 
4780  CLEAR  9  H-2 
4790  RETURN 
5000  REM  * 

5010  REM  * 

5020  REM  *  FUNCTION  FNA 

5030  REM  *  15  MAKE  THE  LABELS  AN 

5040  REM  *  INCREMENT  OF  10 

5050  REM  *  25  AOJUST  LABEL  TO  CORRECT 

5060  REM  *  ENGINEERING  UNIT 

5070  REM  * 

5080  REM  * 

5090  OEF  FNAC85 

5100  J3-0 

5110  FOR  Jl-1  TO  5 
5120  9-B/1000 
5130  IF  B<1  THEN  5155 
5140  J4<J5— J4CJ5-1 
5150  NEXT  J1 
5155  J5-J4CJ5-2 
5157  IF  J5<1  THEN  5185 
5160  FOR  J2-1  TO  J5 
5170 

5180  NEXT  J2 

5185  8-8*1000  9  C-8 

5190  FOR  Jl-1  TO  3 

5200  8-8/10 

5210  IF  B«1  THEN  5240 


5220  J3-J3+1 
5230  NEXT  J1 

5240  B-CIP  <B*10)  "-D  *10*J3 

5245  J-J+l 

5250  FNA-B 

5260  FN  END 

5300  REM  * 

5310  REM  * 

5320  REM  *  TRANSFER  SUBROUTINE 

5330  REM  *  1)  TAKE  IN  TWO  CURVES 

5340  REM  *  AND  MOVE  THEM  TO 

5350  REM  *  REGISTERS  1  &  2 

5360  REM  * 

5370  REM  * 

5380  IF  A>B  THEN  F=B  3  B-A  3  A-F  3  Cl-1  3  C2-2  0  GOTO  53 

90 

5385  Cl -2  0  C2-1 

5390  IF  A-l  THEN  5466 

5420  PCD-PCA) 

5430  FOR  1 1-1  TO  P<D-1 
5440  X(l,  I15-XCA.  ID 
5450  YC1.  I1)-Y(A.  ID 
5460  NEXT  II 
5466  IF  B-2  THEN  5515 
5470  P  <23  — P  (B) 

5480  FOR  1 1-1  TO  P (2) 

5490  X  (2.  ID-XCB.  ID 
5500  Y  (2.  ID  — Y  <B.  ID 
5510  NEXT  II 
5515  GOSUB  6200 

5520  GOSUB  4725 

5530  A-l  3  B-2 

5540  RETURN 
5600  REM  * 

5610  REM  * 

5620  REM  *  PLOT  SUPPLEMENT 

5630  REM  ♦  D  MAKES  V-I  CURVE 

5640  REM  * 

5650  REM  * 

5660  CLEAR 

5670  OISP  "Load  tno  printer  and  press"  3  OISP  "CONTINUE. 
•• 

56e0  PAUSE 

5685  LORG  6  3  CSIZE  <H*  OJ  <23  -U  CD  3  /S83 1*50+ 103 /19.  .  3 
5690  LOCATE  10.  H*  CV  CD -V  <35  3  /8831  *50+10.  10.  H*  <V  CD -V  <33  3 
/883 1*50*10 
5700  FRAME 

5710  SCALE  0.  V1*'/2CC23. 0.  Vl*V2CCl) 

5720  IF  C*-"Y"  THEN  5750 

5730  AXES  V2CC2).  V2CCD.  0.  0.  1.  1.  3 
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5740  GOTO  5760 

5750  GRID  V2CC2).  V2CC1).  0.  0.  1.  1.  1 

5760  FOR  12-1  TO  VI 

5770  HOVE  1 2*V2  CC2)  .  - CV2 CC 1 )  / 1 0) 

5780  LA8EL  VALS CI2*V2 CC2) ) 

5790  NEXT  12 

5800  LORG  8 

5810  FOR  13-1  TO  VI 

5820  MOVE  -CV2CC2)/8).  I3*V2CC1) 

5830  LABEL  VALS  C I P  < 1 3*V2  CC 1 ) -.  55 ) 

5840  NEXT  13 

5850  LORG  6  9  CSIZE  CH*  CV  Cl) -V  C3)  ) /8831*100-»-ia) /16.  .  5 
5860  MOVE  .  S*V1*V2CC2).  -  CV1*V2  CC1) /10) 

5870  LABEL  VOS  CZ1  CC2-1)  +1.  Z1  CC2)  3 »  "  C'i  V4S  CZ  CC2-1 )  +1,  Z  CC 
2) 1 1  ") " 

5880  OEG  9  LOIR  90  9  LCRG  4 

5890  MOVE  -CV1*V2CC2)/10).  .  S*V1*V2CC1) 

5900  LABEL  VOSCZl  CCl-D+l.  Z1  CCl)3j  "  C"»  V4S  CZ  CC1-1) +1.  Z  CC 
1)  3  i  ")  " 

5910  RAO  9  LOIR  0 

5920  N-V1-V2  CC2) / CV  Cl) -V  C3) ) 

5930  0— V1*V2  CC1)  /  CV  Cl)  -V  C3)  ) 

5940  MOVE  CX  C3.  1 )  -V  C3)  )  *N.  CY  C3.  1 )  -V  C3)  )  -0 
5950  FOR  12-2  TO  PC3)-1 

5960  ORAW  CX  C3.  12)  -V  C3)  )  *N.  CY  C3.  12)  -V  C3)  )  *0 
5970  NEXT  12 

5980  CLEAR  9  DISP  "Haw  many  linos  doss  tha  main" 

5990  DISP  "titla  havo?" 

6000  BEEP  200.50 
6010  INPUT  T1 
6020  LORG  4 
6030  FOR  12-1  TO  T1 

6040  DISP  9  DISP  "What  is  lino"iI2j"7"  9  INPUT  T$ 

6050  DISP  9  OISP  "Whora  do  you  wnt  lino"* II  9  DISP  "prin 
tod?" 

6060  DISP  “ CGivo  y-coardinato  from  tha"  9  OISP  "graph.)" 

6070  INPUT  T2 

6080  MOVE  . 5*V1*V2CC2).  T2 

6090  LABEL  TS 

6100  NEXT  12 

6110  GOTO  530 

6200  REM  * 

6210  REM  * 

6220  REM  *  SUPPLEMENTARY  TRANSFER 
6230  REM  *  SUBROUTINE 

6240  REM  * 

6250  REM  * 

6260  IF  XC1.  1)»XC2.  1)  THEN  Rl-2  9  R2-XC1.  1)  9  GOTO  6280 

6270  Rl-1  9  R2-XC2.  1) 


